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1. Introduction 

Recently, there has been a renewed interest in the unusual transport and spectral 
properties of nanotubes, ballistic wires, and quasi-one-dimensional (ID) compounds 
[1, 2]. Quantum effects are strongest at low dimensionality leading to unusual 
phenomena such as charge-spin separation at all energies [2] and persistent currents 
in mesoscopic rings [3]. Thus, the further understanding of the transport of charge in 
low-dimensional correlated systems and materials is a topic of high scientific interest. 

There is numerical evidence of a fundamental difference between the transport 
properties of integrable and nonintegrable ID interacting quantum systems: at finite 
temperatures, T > 0, the integrable systems behave as ideal conductors in the metallic 
quantum phases and as ideal insulators in the insulating phases, with the concepts of 
an ideal insulator and conductor defined in Ref. [4]. In contrast, the nonintegrable 
ID interacting systems are generic conductors and activated ones in the metallic and 
insulating phases, respectively. While in the trivial case of ID integrable systems whose 
Hamiltonians commute with the current operator the ideal insulating and conducting 
behaviors are easy to confirm, there is the expectation that such ideal behaviors might 
also occur in ID integrable quantum systems whose Hamiltonian does not commute with 
that operator, such as the ID Hubbard model [5, 6, 7]. However, the studies of Ref. 
[8] rely on the generalization of the thermodynamic Bethe-ansatz equations introduced 
in Ref. [7] to the model in a presence of a vector potential [9] and seem to reveal that 
for half filling the ID Hubbard model does not remain an ideal insulator for T > 0, in 
contrast to the general predictions and expectations of Ref. [4]. 

Solvable lattice models such as the ID Hubbard model [5, 6, 7, 10] and the XXZ 
chain [11] are often used as toy effective models for the study of the unusual properties 
of quasi-lD compounds [12, 13]. Although the ID Hubbard model was diagonalized 
long ago [6, 7] by means of the coordinate Bethe ansatz (BA) [14, 15], the involved form 
of BA wave functions has prevented the full calculation of dynamic response functions. 
The study of the asymptotic of correlation functions and of the low-energy dynamical 
properties was performed by combining the BA solution with other methods, such as 
conformal-field theory [16], bosonization [17, 18], the pseudo-particle formalism[19], and 
scaling methods [20]. 

In this paper, we study the charge and spin currents carried by the elementary low- 
energy and finite-energy excitations of the ID Hubbard model. The paper is organized 
as follows: In Sec. II we summarize the basic information about the model and the 
BA solution needed for our study; The energy and current spectra of the elementary 
excitations is the subject of Sec. III. Finally, in Sec. IV we present the concluding 
remarks. 
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2. The model and its Bethe-ansatz solution 

The Hamiltonian of the ID Hubbard model on a periodic iV a -site chain reads, 

H = "* 5>5,*9*+i,* + c} +1>ff c,- a ) + Uj2(n^ - l/2)(n u - 1/2) , (1) 

3i u 3 

where the operator c] (and Cj a ) creates (and annihilates) an electron of spin projection 

' t 

a at the site of index j = 1, 2, 3, N a and n^ c = cj a Cj :(T is the number operator at the 
same site. We use units of lattice constant one such that L = N a , where L is the system 
length. We denote the electron number and the spin-projection a electron number by 
iV and N a , respectively, such that N = [iVj + iVjJ. Moreover, we denote the states spin 
and 77-spin values by S and rj, respectively. 

The model as written in Eq. (1) has both a spin and a rj-spin SU(2) symmetry 
[21, 22, 23, 24]. The generators of the 77-spin symmetry are given by, 

L L 

3=1 3=1 
1 L 1 

V z = ^J2^i + n ^ ~ 2 1, 

3=1 

[77, rf] = -2r] z , [77, r] z ] = rj, [r]\r] z ] = -rf . (2) 

The global symmetry of the model (1) corresponding to these two SU(2) symmetries is 
5*0(4), since half of the irreducible representations of SU (2) (^) SU (2) are excluded. The 
BA solution refers to the Hilbert subspace spanned by the lowest-weight states (LWSs) 
of both the spin and r^-spin algebras. The subspace which is not associated with such a 
solution is spanned by the energy eigenstates obtained by applying onto the LWSs one 
of the off-diagonal generators of the corresponding two algebras [25]. 
The charge and spin current operators of the model read [19], 

L 

Jp = -eit 5Z( c i- c i+!- ~ c i+i<x c i-) 

a j = l 
1 L 

■J az = Yl Yl (J (4- C i+ 1 - - c ]+ia c i-) • ( 3 ) 

<T j = l 

In order to calculate the expectation values of the charge and spin current operators 
it is convenient to consider a uniform vector potential A x e x , which modifies the hopping 
term along the chain by the usual Peierls phase factor, t — > texp(±i(p a /L). Following 
such a procedure, the Hamiltonian becomes, 

3,<r 

+ ^K r l/2)K r l/2). (4) 

3 
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For a given energy eigenstate \m), the charge and spin current expectation values 
Jm = ('m\J p \'m) and = (m\J az \m), respectively, can be expressed as follows [26], 



dE m ((f)) 



d{<t>/L) 
dE m (4>) 



= 0T 

0=0 



J ™ d(<P/L) 



= 0T = ~<f>l- ( 5 ) 

=0 

The (pa > Hamiltonian (4) remains integrable and can be diagonalized by means 
of coordinate BA [9]. One can introduce two generalized SU(2) symmetries for the 
(f) a > case [27]. Since 2i] and S remain good quantum numbers, one finds that the 
BA solution refers to the LWSs of both the 77-spin and spin generalized algebras. Thus, 
all the energies of the tower of states such that 2r\ > [L — N] and 2S > [JV-p — N{\ have 
the same energy as the corresponding LWSs. It follows that the BA numbers can be 
related to the values r] and S of the states of each 7/-spin and spin tower, respectively. 
Since the studies of the ensuing section refers to both LWSs and non-LWSs, here we 
express the sum rules of the BA numbers in terms of the good quantum numbers r\ and 
S. Moreover, we provide the simplified expressions in terms of the electronic numbers 
N and N a which correspond to the LWS of each tower only. 

The solution of the Hamiltonian (4) by the BA leads to the following equations 
[9, 28], 

M . , . 

e ikjL = TT sin kj — Ap + m 
^■jL sin kj — A/3 — iu ' 

e ^i-0t)ff A^sinfcj+m 

J- J- A^ — sin kj — iu 

j = i 1 j 

n l A 7 - A^ + i2u 
=i A^-A f3 -i2u- { ] 

Here and throughout this paper u = U /At, the numbers N c and M such that < N c < N 
and < M < N[, respectively, are defined below, and A is the spin rapidity [7]. 
Takahashi string hypothesis states that in addition to the real solutions for A 7 and kj, 
there are solutions involving complex kj and A 7 values. The spin string As of legth n is 
characterized by [7], 

A? = A™ + (n + 1 - 2j)iu. j = 1, 2, . . . , n. , (7) 

where A™ is the real part of the complex number. The charge k — A string of length n 
includes 2n ks and n As such that, 



A^ = A' n y + (n + 1 - 2j)iu, n = 1,2, . . . ,n, 
k^ = n - sin -1 (A'™ + niu), 
P = sin" 1 (A'! + (n-2)m), 
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kt, = sin" 1 (A'; + (ra-4)m), 

kry = 7T k -,, 



. . . , 

A;f- 2 = sin- 1 (A / ;-(n-2)2 M ), 

A; 2 ™ = 7 r-sin- 1 (A / "-mu). (8) 

By use of Eqs. (7) and (8) in Eqs. (6) we arrive to the following transcendental equations 
[28, 29], 

00 Mn eink- — A™ °° M ' n dn^-A" 1 

^ = an, + * - £ £ Oi^y^) - E E "f^ 1 ^ 1 ) 

n=l /3=1 n=l /3=1 

L(sin _1 (A / ™ + inu) + sin -1 (A'™ — inu)) = 2irJ'™ — + 

N—2M' \>n _ ■ i A /n - A /m 

+ E *( 7 V E e ^-^ 

rat u 

j=l m,/3 

x 2U ' , A™ — sin fc,- ^ A" - AT 



E ^ 7 - ' ) = 2 ^7 + - ft) + E ©n m (-^ L ) , (9) 

where 9{x) = — 2tan(x) and 



nu ' ^-^ u 



e nm (x) = 9 I )+29( , X —— )+■■■ + 29 ' 



n — m\J \\n — m\+2J \n + m — 2 

+ 6 ( — - — ) , for : n ^ m 

\n + m J 



29 : 

\2, 



;) + »o + - + »(5^) + '(s)- 

for : n = m . (10) 

Here {/-,-, J'™, J™} are the actual quantum numbers whose values define the energy 
eigenstates and thus determine the energy and current spectra of the elementary 
excitations studied in the ensuing section. Following the notation of Takahashi [7], 
we introduce the numbers, 

M c = E nM n> M s = Y,nM n . (11) 

n n 

In these expressions M n and M' n are the numbers of spin A strings of length n and charge 
k — A strings of length n, respectively. The values of N c and M are then uniquely defined 
by the following sum rules, 

N c = L - 2 V - 2M C , M = M c + M s = ^-r]-S, (12) 

which for a LWS of the 77-spin and algebras such that 2r\= [L — N] and 2S = [N^ — N±] , 
respectively, simplify to, 

N c = N — 2.1/,. M = N i . (13) 
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The above quantum numbers {Ij, J'™, J™} can be integers or half-odd integers (HOI) 
according to the following prescriptions: Ij is integer (or HOI) if ^2 m (M m + M^) is even 
(odd); J™ is integer (HOI) if N c -M n is odd (even); J'™ is integer (HOI) ilL-(N c -M' n ) 
is odd (even). They have values in the following ranges, 

T L 

h < j ' 



J'; < -(L - TV + 2M C - J2 tnmM'J, 

m=l 

1 °° 

J 7 " < -(TV - 2M C - ^ t nm M m ) , (14) 



m=l 



where t nm = 2 min(n, to) — 8 nm . 



The energy and momentum spectra are given in terms of the BA quantum numbers 
as follows, 

E= -J2 2tcosk i + AtlZ\J 1 — (A'™ — inu) 2 



j'=i 



P = 



|(E^E^) + E(-f^) 



1 



+ 7r (M c + 77--[L-iV]), (15) 

where 1Z refers to the real part and rc(M c + rj — ^[L — Nj) simplifies to irM c for a LWS 
of the 77-spin algebra. 

All energy eigenstates associated with the BA solution are described by different 
occupancy configurations of the quantum numbers appearing in the system of coupled 
equations given in Eq. (9). For example, for the ground state there is no complex 
solution for Eqs. (6) and Ij and J 7 are successive numbers centered around zero. Thus, 
the quantum number occupancy configuration for even N c = N = L and odd N c /2 
corresponds to, 

TV — 1 TV-3 TV — 1 
j ~ 2~' 2~' " ' ' ' 

M-l M-3 M — 1 . . 

J, = 2->-~,—2-. (16) 

In some of the figures presented in the ensuing section we measure the energy 
relative to the ground-state energy. Such a choice corresponds to the following general 
energy spectrum, 

— inu) 2 

j=l n,a 

- I [TV C + 2M C - 1] - 2/i(TV a — TV) — 2/i # (TV T - TVJ , (17) 
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where \i is the chemical potential, no the Bohr magneton, and H the magnetic field. 

Although the Bethe-ansatz equations (9) refer to L » 1, we have used these 
equations in the numerical study of finite-L chains and obtained results for several 
quantities in excellent numerical agreement with the known exact values. Thus, in 
the ensuing section we use these equations to derive numerically the charge and spin 
currents carried by the elementary excitations for finite values of L, as well as their 
energy spectra. The obtained results are a good approximation for the corresponding 
current spectra of the L » 1 quantum problem which the equations (9) refer to. 

3. Current spectra for a finite-size system 

3.1. The case of half-filling 

The zero-magnetization and half-filling ground state is both a spin and r^-spin singlet. 
Thus, the simplest elementary excitations are spin-triplet excitations such that (77 = 
0, S — 1), spin-singlet excitations such that (77 = 0, S — 0) whose energy spectrum is 
degenerated to that of the spin-triplet excitations, charge 77-spin-triplet excitations such 
that (77 = 1,S = 0), and charge ?7-spin-singlet excitations such that (77 = 0, S — 0) 
whose energy spectrum is degenerated to that of the charge 77-spin-triplet excitations. 
The energy spectra of these four elementary excitations have been previously studied 
[30]. For a comparison, we evaluate the energy spectrum and the charge and spin 
currents of all the energy eigenstates corresponding to the above four types of half- 
filling elementary excitations. The energies considered in this subsection correspond to 
the general energy spectrum provided in Eq. (15). 

The half-filling and zero-magnetization ground state considered in this subsection 
is characterized by charge and spin distributions given by, 



where J is the Bessel function of zero order. Following standard BA procedures, the 
evaluation of the energy, charge-current, and spin-current spectra studied below involves 
the use of corresponding distributions for the elementary excitations. 

Spin-triplet excitations. Such elementary excitations are obtained by introducing 
"holes" in the spin distribution of the numbers J 7 relative to the ground-state occupancy 
configuration. For each value of the excitation momentum and energy there is a spin 
tower of three S = 1 states, differing in the spin projections 0, ±1, but all having the 
same values, 




(18) 



N C = N = L, M = L/2-l, 



(19) 



and 



N-2 



N -4 



N 
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M + l Af-3 M + l 
^- 2 , 2 2 ' 1 j 

for the BA numbers. Hence, there are two holes A^ and A 2 in the spin distribution. The 
BA equations become, 

kjL = 2nlj + T - 2 2^ tan" 1 3 - '- 



u 

0=1 



N \ ■ 1 

Ei — sm kj j 
2 tan" 1 ^ 3 - = 2vr J 7 + - T ) 



+ 2 J] tan" 1 (21) 

/3=1 

Following the usual procedures of BA [7, 28], one then introduces the following charge 
and spin distributions, 



p{k) — — + cos/c j ^(sin k — X)a(X)dX, 
2tc J 

cr(A) + a h {\) = I K ± (X - sin k)p{k)dk 



-Jk 2 (X- X')a(X')dX', (22) 

where a h (X) = [5(X — A^) + 5(A — A^)]/^ and K n (x) = nu/[ix{n 2 u 2 + x 2 )}. In the presence 
of the flux the momentum deviation corresponding to this elementary excitation is given 
by, 

AkjL = 0| 
L/2-1 

— 2n ^(sin kj — Xg)(cos kjAkj — AXp), 

0=1 

AT 

27T i^i(A 7 — sin fcj)(AA 7 — cos kjAkj) = — 0|) 
j'=i 

L/2-1 

+ 2tt ^ K 2 (A 7 - A^)(AA 7 - AA^) . (23) 
0=1 

Use of Eq. (22) then yields, 

Akp(k) = -p- + I K^sink- X)AXa(X)dX, 
2ixL J 

AX[a{X) + a h {X)\ = _ f _ A')AAV(A / )dA / 

2ixL J 

+ J Xi(A- sin k) cos kAkp{k)dk . (24) 
The corresponding energy deviation is given by 

AE((f))=2tL J sin kp(k)Akdk . (25) 




Figure 1. The half-filling energy spectrum (left) and spin current spectrum (right) 
of the spin triplet excitations for u = 10, TV = L = 66, and M = 33. 



Our next task is the solution of the equations given in (24). Inserting the result obtained 
for Akp(k) into Eq. (25) we find, 

"sech [^(sin/c - Af)] sech [^(sinA; - A£)] 



AE 



t(4>i - <h) 

8nuL 



dk sin k 



*o(Af) 



a (\ h 2 ) 



where we have used the relation AA ~ (0^ — 0j)/(47rLer o (A), which was obtained by- 
Fourier transformation. It then follows that when <pi = 0j the charge current defined 
by Eq. (5) vanishes. 

In general case, X h > 1, Eq. (26) can be solved approximately with the result, 



r 



2u 



tanh 



~2u 



+ tanh 



7tA| 

2u 



(27) 



Since the momentum carried by the spin elementary excitation has the form q = 
2 tan -1 e~ A / 2 " — | and the quantum number of the charge part changes from half-integer 
to integer, what generates a momentum shift n, the spin-current spectrum has the 
following form, as was also observed in the antiferromagnetic Heisenberg model [31], 

J = — [smq 1 + smq 2 \, 

q =n + q 1 + q 2 , q 1 , q 2 G [- J, J] . (28) 



Here q\ and q 2 are the momenta of the two spin-distribution holes and q denotes the 
total excitation momentum. In figures 1 and 2 we plot the energy and spin-current 
spectra of the spin-triplet excitations for u — 10 and u — 1, respectively. Such spectra 
were obtained by solving numerically the BA equations. 



(26) 



0.5 



1 

q/n 



1.5 




Figure 2. The half-filling energy spectrum (left) and spin current spectrum (right) 
of the spin triplet excitations for u = 1, N = L = 66, and M = 33. 



The group velocity v a ^ (q) and the effective spin (charge) e a ^ are defined as, 

In this equation a and p denote the spin current and charge current, respectively, and 
e = — 1 and e = 1/2 for the charge and spin cases. Thus, the group velocity of a single 
spin-distribution hole reads, 

tix 

v°{q) = -—smq, (30) 
whereas the corresponding effective spin is given by, 

** = ~. (3D 

Note that the total spin current is positive. 

Spin- singlet excitations. The second type of elementary spin excitations studied 
here corresponds to the spin-singlet excitations whose energy spectrum is degenerated 
with that of the spin-triplet excitations considered above. Such excitations have a spin 
string of length one. Thus, the BA numbers are given by, 

M 1 = M-2, M 2 = 1 , (32) 



and 



t2 



N -2 iV-4 N 



2 ' 2 ' ' 2 ' 
M - 1 M - 3 M - 1 



4 = . (33) 
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Figure 3. The half-filling energy spectrum (left) and spin current spectrum (right) 
of the spin-singlet excitations for u = 10, N = L = 66, and M — 33. 
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Figure 4. The half-filling energy spectrum (left) and spin current spectrum (right) 
of the spin-singlet excitations for u = 1, N = L = 66, and M = 33. 



It follows that there are again two holes in the spin distribution, A^ and A£. The BA 
equations are given by, 

L/2-2 . 

Sin Kj — A/3 



kjL = 2-kIj + T - 2 tan" 

/J=i 
_ x sin fcj — A 



u 



— 2 tan 



2m 




Figure 5. Scaling analysis of the maximum current carried by spin-singlct excitations 
for u = 1 (left) and u = 10 (right) respectively. 



N 



2 1^- 1 A7 = 2tt J 7 + (0, - T ) 



L/2-2 



tan 



-1 A 7 — A/3 

2m 



+ 2 tan" 1 + 2 tan" 1 



3m 



N 

2^ tan 

3=1 

L/2-2 

P 



'^i = 2,j! 2 >(=0) + 2to^ l) 



i A - A/3 i A - Afl 

tan -1 - + tan" 1 p 



u 



3u 



(34) 



where A denotes the rapidity of the spin string excitation of length two. The deviations 
of the charge and spin distributions are such that, 



Akp(k) = + / K 1 (sink - \)A\a(\)d\ 
2%L J 



+ -iG (sin A; - A)AA 



AA[cx(A) + a h (X)] 



01 - 0T 
2ttL 



- J K 2 (\ - X')AX'a(X')dX' 

+ J cos k/\kp{k)Ki{\ — sin k)dk 

-i[X 1 (A-A) + K 3 (A-A)]AA 
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AA<r 2 (A) 



- J [Ki(A - A) + K 3 (A - X)]AXa(X)dX 
+ J cos kK 2 (A-smk)Akp(k)dk 

- J [2K 2 (A - A') + K A (A - A')]AA<7 2 dA . (35) 

Here a 2 (A') = 5(A' - A)/L and a h {\) = [5(X - Af) + 5(X - Af)]/L. Integrations in the 
variables fc, A, A lead to, 



cos kAkp(k) = 
y AAa(A)rfA + -^AA = ^ ^^I dX- J AXa h (X)dX, 

jAA + i y AAcr(A)dA = J \~^ dX . (36) 

Finally, the deviation associated with the spin string excitation of length two has the 
form, 



AA = (AAf + AA^) + J Vi 4n dX , (37) 

where the second term on the right-hand side can be omitted because it does not 
contribute to the spin current. Use of this expression in the first equation of (35) 
leads to, 

Akp(k)^^-+ [ K^sink- X)AX[a + a h ]dX. (38) 
2nL J 

Since a (X) a(X) + a h (X), by means of the same procedure already used for spin-triplet 
excitation, we find that to first order the energy deviation induced by the external flux 
vanishes. Thus, both the charge and spin currents carried by this type of elementary 
excitation vanish. 

In figures 3 and 4 we plot the energy and spin current spectra of the spin-singlet 
excitations for u — 10 and u — 1, respectively. The two small features of the spin- 
current spectrum result from finite-size effects and disapear in the thermodynamic limit, 
as shown in figure 5. 

We emphasize that although the group velocity of the two spin-distribution holes 
of the spin-singlet excitations is finite, the corresponding effective spin vanishes. 

Charge rj- spin-triplet excitations. For each value of the excitation momentum and 
energy there are three types of such 77 = 1 elementary excitations, which correspond to 
the three values 0, ±1 for the 77-spin projection. All these excitations have again the 
same BA numbers, 

N = L - 2, M = N/2, (39) 



2-2 2 

ql% q/K 

Figure 6. The half-filling energy spectrum (left) and charge-current spectrum 
(right) of the charge 77-spin triplet excitation for u — 1 0, N = L = 46, and M = 23. 




and 

I=- L -... L --X 

3 2 ' '2 

M-l M- 1 
J, = — 2~ • ( 4 °) 

There are two holes A;^ an d ^2 i n the charge distribution such that, 

p h (k) = j[5(k-k h l ) + 5(k-k h 2 )}. (41) 
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Thus, the distributions p{k) and er(A) satisfy the following equations, 

p(k) +p h {k) = ^- + cosk [ K^sink - X)a(X)dX 
2n J 

cr(A) = J K 1 (A-sinA;)p(A;) 

- J K 2 (X - X')a(X')dX' . (42) 

The energy deviation corresponding to the two charge-distribution holes read, 

AE((f)) = -2t [sin k^Ak^ + sin k%Ah%] , (43) 
where the momentum deviation Ak is approximately given by, 

A A; « . (44) 

2nLp (k) v y 

It follows that the charge current spectrum is such that, 

sin k\ sin k 2 



J p = -2t 



Lpo(^) p (k h 



2, 



(45) 



while for k h = —k it reads, 

sin k] sin A> 

J p = 2t 



_Po(h) Po(k 2 

k =k 1 + k 2 . (46) 

In turn, the spin current vanishes. In the strong coupling limit, u = U/4t » 1, one 
has that po(k) ~ l/2n and thus the spectrum simplifies to, 

J p oc sin ki + sin k 2 . (47) 

In this case the corresponding energy spectrum can be expressed as the sum of three 
cosine functions, in addition to the energy gap. Hence, the velocity of a single charge- 
distribution hole simplifies to v p = —2t sin kh, whereas the effective charge is — e. This 
corresponds to a positive current in units of e. 

Again, we used the BA equations to calculate the energy, charge-current, and spin- 
current spectra for u — 10 and u—1. (The spin current vanishes for the charge ^-spin- 
triplet states considered here.) The energy and charge-current spectra are plotted in 
Figs. 6 and 7 for u = 10 and u — 1, respectively. Note that the charge-current spectra 
features have a stronger linear character for u — 1 than for u — 10. We interpret this 
effect as due to the weak-coupling peak in the charge distribution as a function of the 
momentum k. When a hole is created away from zero momentum it is less affected by 
the other charges, leading to J p oc q. 

Charge rj- spin- singlet excitations. These r\ = elementary excitations are those 
whose energy spectrum is degenerated with that of the charge ^-spin-triplet excitations 
considered above. Such 77-spin-singlet excitations contain one charge string of length 
one. The BA numbers are then given by, 
N 

M 1 = --l, M[ = l, (48) 
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Figure 8. The half- filling energy spectra of charge 77-spin singlet excitation for the 
strong coupling u = 1 (left), and weak coupling u = 10 (right). 
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J( = . (49) 

These excitations involve two holes k\ and k\ in the charge distribution. The BA 
equations are such that, 

L/2-1 . 
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(50) 



where A is the rapidity involved in the charge string of length-one. Moreover, we find, 



Akjip + p 11 ) = 1 P-+ [ ATi(sinfc- A)AAa(A)dA 
27rL J 




+ -K^smk - A)AA, 



AAcr(A) 
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2ttL 

J K\(h — sin fc) cos kAkp{k)dk 
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(51) 



where o"'(A') = <5(A' — k)/L. Integrations involving the variables A, k, and A lead to. 



This is consistent with p h = [5{k - fcf) + 5{k - k%)]/L and A = (sinA;^ + sin£;£)/2. 
The energy of the present elementary excitation is, 



Thus, the corresponding energy deviation in the presence of the external flux can be 
expressed as, 



Note that the momentum of the two charge-distribution holes is determined by the value 
of the rapidity A. The important point is that their contributions to the charge-current 
cancel in the thermodynamical limit. 

Numerical solution of the BA equations with the quantum- number occupancy 
configurations given in Eq. (49) leads to the energy spectra plotted in Fig. 8. From use 
of Eq. (5) we could confirm that both the spin and charge carried by these elementary 
excitations vanish in the thermodynamical limit [27]. 

3.2. Half- filling with non-zero magnetization 

While the energy spectra plotted in figures 1-8 correspond to the general energy 
spectrum given in Eq. (15), in the figures 9-14 presented below the energy spectra 
of the elementary excitations refer to the general spectrum provided in Eq. (17). 

Here we consider that the initial half-filling ground state corresponds to a finite 
spin density. That is achieved by the presence of a magnetic field. Thus, for such a 
ground state and corresponding elementary excitations one has that M < N/2 and the 
integration limit of the spin variable A changes from oo to a finite cut-off. For odd values 
of M, the distribution of the corresponding spin quantum numbers is still of the form 
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Figure 9. The half-filling energy spectrum (left) and spin-current spectrum (right) 
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Figure 10. The half-filling energy spectrum (left) and spin-current spectrum (right) 
of the spin excitation with a string of length two for u = 1, N = L = 75, and M = 25. 



given in Eq. (16). The charge excitations do not change much due to the spin-charge 
separation. For the spin sector, however, the excitation spectrum presents qualitative 
changes. 

There are still excitations which increase the value of S by one. However, such 
excitations require a minimum finite energy and thus the corresponding excitation 
energy spectrum has a gap. The spin distribution of these elementary excitations 
displays two holes. Moreover, such elementary excitations carry in general a finite spin 
current. Additionally, since the integration limit of A is now finite, the external flux can 



shift the whole integration region. That process also contributes to the spin current. The 
corresponding current feature is expected to be a linear function of the momentum q. 
The same kind of spin-current feature arises now from elementary excitations involving 
one spin string of length two. For the case of the zero-magnetization ground state such 
excitation is of spin singlet character. However, in contrast to the zero-spin case now the 
two spin-distribution holes contribute to the spin current. Consideration of elementary- 
excitations whose deviations from the quantum number configuration occupancies is 
the same as that of equations (33) for the spin singlet excitation, leads in the case of 
an initial ground state with finite spin density to the energy and spin-current spectra 
plotted in Figs. 9 and 10 for u = 10 and u — 1, respectively. Note that now such 
excitations lead indeed to a finite value for the spin current. 



3.3. Away from half filling 

Finally, let us consider that the initial ground state is metallic and thus refers to an 
electronic density away from half-filling. The ground-state spin density is considered to 
be zero. In this case, besides the ground state, there also exist real rapidity solutions 
for some of the elementary excitations. For simplicity, we still consider the case of N/2 
odd but now with N < L. The simplest elementary excitation away from half-filling 
corresponds to removing one Ij from the ground-state occupied charge-distribution 
sector and adding a new one outside such a sector. Such an excitation has a "particle- 
hole" character. It is characterized by the following BA numbers, 

r n-i N-i 

+ 71 + 1,...,—— ,I n }. (55) 



2 7 7 2 

where \I n \ > (N — l)/2. We consider excitations such that the A spin distribution 
remains unchanged. 

In figures 11 and 12 we plot the energy and charge- current spectra of such 
elementary excitations for u — 10 and u — 1, respectively. Note that for strong coupling 
(u — 10) the charge current may have negative values whereas for weak coupling (u — 1) 
it has always positive values. 

In the strong coupling limit the charge BA equation simplifies to, 

kjL = 2-Klj + T . (56) 
It follows that the energy deviation can be written as, 

A£(0) = 2t sin kjAkj . (57) 

3 

Furthermore, the current reads, 

J p = 2t{sm(2k F + q)- sin(2A; F )] , (58) 

where 2kp = nN/L denotes the "Fermi momentum" and q the momentum. Note that 
Eq. (58) provides the charge current for an excitation where a charge is removed from the 




Figure 11. The quarter-filling energy spectrum (left) and charge current spectrum 
(right) of a charge "particle-hole" excitation for u = 10, L = 100, TV = 50, and 
M = 25. 




Figure 12. The quarter-filling energy spectrum (left) and charge current spectrum 
(right) of a charge "particle-hole" excitation for u = 1, L = 100, N — 50, and M = 25. 



"Fermi point" and created outside the ground-state "Fermi sea". As long as 2k p > tt/2, 
the current in the vicinity of 2k f is negative. On the other hand, the effective charge 
carried by this excitation is negative. 

The contribution to the charge current of the charge-distribution hole excitation is 
straightforward to obtain and reads, 

J p h = 2t[sin2k F - sin k h ], k h G [-2k F , 2k F ] , (59) 



where kh is the momentum of the charge-distribution hole. At 2k p — ir/2 it leads to a 
feature in the charge current spectrum defined by the function [1 — cosg/2]. Combining 
the "particle" and "hole" contributions one finds, 

JP = 2t[sin{2k F + k) - sin k h ] . (60) 

The group velocities of the charge "particle" and "hole" are such that, 

v p v = —2t sin hp, v p h = 2t sin k h . (61) 

Hence, the effective charge carried by the "particle" and the "hole" are e(= —1) and 
— e(= 1), respectively. We recall that all this analysis applies to the strong coupling 
limit only. 

In turn, in the weak coupling limit the interaction between the charge holes is 
so weak in the "Fermi sea" that the charge current is almost a linear function of the 
momentum. For the "particle" at quarter filling, however, since the distribution function 
p(k) becomes a very narrow peak and the "Fermi surface" is compressed, the charge 
current changes from negative to positive. 

We have shown above that both the spin and charge currents carried by the charge 
r7-spin-singlet excitations containing one charge string of length one vanish. However, 
away from half-filling the corresponding charge excitations containing one charge string 
of length one have different properties. The quantum number occupancy configuration 
of such elementary excitations is still described by Eq. (49). We have solved the 
corresponding BA equations. The found energy spectrum and charge-current spectrum 
are plotted in Figs. 13 and 14 for u — 10 and u — 1, respectively. In contrast to the 
corresponding excitations relative to the half-filling ground state, the present excitations 
carry charge but no spin. 

4. Concluding remarks 

In this paper we have studied the spin and charge currents carried by the elementary 
excitations of the ID Hubbard model. Most of our results refer to half-filling. Both the 
charge ^-spin-singlet and spin-singlet elementary excitations considered in our study 
carry no charge and no spin. Moreover, the spin-triplet excitations considered in this 
paper carry spin but no charge, while the charge ^-spin-triplet elementary excitations 
carry charge but no spin. 

Our results reveal that the present quantum liquid is not an ideal insulator for half- 
filling, in apparent contradiction to the general predictions of Ref. [4]. Indeed, ideal 
insulating behavior requires that all half-filling states carry no charge current. However, 
we note that out of the three types of charge 77-spin-triplet elementary excitations of 
77-spin projections 0, ±1, the states with projection are indeed half-filling states. Our 
results show that such charge 77-spin-triplet states carry finite charge current both in 
the strong and weak coupling limits. According to the results of Ref. [4], this implies 
a finite value for the charge stiffness D(T) at finite temperatures T > 0. This result 
agrees qualitatively with the studies of Ref. [8], which lead to a finite value for D{T) at 




Figure 13. The quarter-filling energy spectrum (left) and charge-current spectrum 
(right) of a charge excitation containing one charge string of length one at for u = 10, 
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Figure 14. The quarter-filling energy spectrum (left) and charge-current spectrum 
(right) of a charge excitation containing one charge string of length one at for u = 1, 
L = 100, N = 50, and M = 25. 

half filling and T > 0. However, such studies did not take into account the r\ > and 
i] z = states which carry the charge current at half filling and lead to D(T) > for 
T > 0. Moreover, both the studies of the present paper and the related studies of Ref. 
[27] provide strong evidence that the rj = and i] z = states considered in Ref. [8] do 
not carry charge current and do not contribute to D(T). 

The model studied here is integrable and according to the general arguments of Ref. 
[4] should display ideal insulating behavior at half filling. That issue is clarified elsewhere 



by study of the microscopic mechanism which is behind the half-filling properties 
concerning transport of charge [27]. 

Finally, our results reveal the occurrence of charge-spin separation at finite energies, 
since some of the elementary excitations studied here have an energy gap relative to 
the initial ground state. Such a finite-energy charge-spin separation deserves further 
studies. We note that the photoemission studies of Ref. [2] have detected a spin-charge 
separation in quasi- ID compounds for the whole energy bandwidth. 
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